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ABSTRACT 

Recent simulations of merging black holes with spin give recoil velocities from gravitational radiation 
up to several thousand km s^ 1 . A recoiling supermassive black hole can retain the inner part of its 
accretion disk, providing fuel for a continuing QSO phase lasting millions of years as the hole moves 
away from the galactic nucleus. One possible observational manifestation of a recoiling accretion disk 
is in QSO emission lines shifted in velocity from the host galaxy. We have examined QSOs from the 
Sloan Digital Sky Survey with broad emission lines substantially shifted relative to the narrow lines. 
We find no convincing evidence for recoiling black holes carrying accretion disks. We place an upper 
limit on the incidence of recoiling black holes in QSOs of 4% for kicks greater than 500 km s _1 and 
0.35% for kicks greater than 1000 km s _1 line-of-sight velocity. 
Subject headings: galaxies: active — quasars: general — black hole physics 



1. INTRODUCTION 

Recent breakthroughs in numerical relativity have 
allowed the possibility of simulating in complete general 
relativity the final orbits of a binary black hole up to 
and past the merger. Analytical approximations have 
for some time indicated that the merger of unequal 
mass black holes will, through anisotropic emission of 
gravitational radiation, impart to the final black hole a 
substantial 'kick' of up to mor e than 1000 kilometers 
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per se cond (iFitchett 
120061 ISopuerta et"alT 
and references therein). Numerical simulations have 
only begun to explore the parameter space of black 
hole mass ratio, spin magn itude, and spi n -orbit 
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orientations in such mergers ([Herrmann et alj 120061: 
Bak er et alJl2006t IGonzalez et al.ll200"7T l Herrmann et al 



are 
l 



20071 iKoppitz et all 120071 : iCampanelli et al 
2007t IGonzalez et all l2007t iBaker et all l2007t 
Tichv fc Marronettil l2007h . Interestingly, they 
consistently showing kicks of order 100 to 1000 km s 
depending on black hole spin inclinations, up to a 
maximum of ~ 2,500 km s _1 for spins anti-aligned 
and perpendicular to the orbital angular mom entum 
(jGonzalez et all 120071; [Tichv fc Marronettil I2007D . Such 
large recoil velocities have significant astrophysical 
implications for galactic mergers, since even velocities 
of order 1000 km s _1 can be greater than the escape 
velocity of moderate-sized elliptical galaxies and spiral 
bulges, and much greater than the <300 km s" 1 escape 
velocity for dwar f galaxies (jCampanelli et al.l 120071 : 
iMerritt et al.|[2004l and references therein). 

Super-massive black holes (^ 10 8 M(T)) will be forme d 
during the merger of galaxies (Begelm an et al.l Il980f ). 
The binary orbit will decay quickly as a result dynam- 
ical friction due to the stellar background. The orbit 
may stall at a radius ~ 1 pc, but thi s may be overcome 
by the presence of a nuclear gas disk (|Escala et alJl2004t 
iDotti et al.l 120061 : ISesana et al] 120071) . For a black hole 
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merger taking place in an AGN, the accretion disk will 
remain bound to the recoiling black hole inside the radius 
R = 1.33 x IO^Ms/i^qqq cm where the orbital velocity 
is equal to the recoil velocity. Here, Mg = M/10 8 M© 
and wiooo = w/1000 km s _1 . The retained disk mass 
in such a case, assumin g an a disk (|Shakura fc Sunvaevl 
[197! iFrank etlifll200"2T ) will be: 



M(R) = (10 bbu M Q )a 



- 4/5 M s 1/4 M n 7/10 i?^ 4 / 14/5 



or 



M(v) = (10 8 - 02 M Q )a 



4/5 M s 3/2 M n 7/lt \ 



-5/2 
1000 



(1) 



(2) 



where Mq is the accretion rate in sol ar ma s ses pe r year. 
Stability requires M disk < M B h (see lLoebl (|2007l )). 

For a black hole binary contained in an accretion 
disk, the two holes will empty out a 'gap' with a 
radiu s of approximately twice the b inary semi-major 
axis (|Macfadven fc Milosavlieyidl2006l) . This will refill 
quickly after the kick ( Loebl l2007h . and QSO activity 
will resume. The disk mass will be sufficient to fuel QSO 
activity over a disk consumption time w M(v)/Mq rs 



(10 8 yr)alf /5 M fi 3/2 M n 
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" J 1000 



. This 'wandering QSO' 



phase could last for a time comparable to the pre-merger 
phase, and would result in either a QSO displaced a num- 
ber of kpc from the galactic nucleus or QSO emission 
lines shifted relative to galactic systemic velocity. 

Obse rvations of nearby AG N do not show displaced 
nuclei (|Libeskind et alj .2006). This may be one indica- 
tion that large kicks rarely occur during an active AGN 
phase. 3 Alternatively, these kicks may be observed in the 
velocity of the AGN emission lines. The broad emission- 
line region (BLR) of QSOs corresponds to radii within 
which the disk will remain bound to the post-merger 
black hole. The BLR dynamical timescale is ~ 10 2 years, 
so the BLR should be regenerated quickly from the disk 
following the recoil and resumption of AGN activity. The 
narrow emission lines, in contrast, arise from gas pre- 
dominantly orbiting in the potential of the host galaxy 



3 The quasar HE0450-2958 ( Mag ain et"all [2005) has been 
suggested as such a candidate, but questi ons remain about 
whether an ejected black hole is ind icated ( Hoffman & LocB l2006l ; 
IMerritt et a l.' 2006; Kim et al. 20071). 
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that will not follow the recoiling black hole (|Merritt et al.l 
120061) . The displaced QSO will still ionize the interstellar 
gas, producing narrow emission lines, albeit different in 
detail from a normal narrow line region (NLR). There- 
fore, the broad emission lines associated with a recoiling 
disk will appear shifted with respect to the galaxy sys- 
temic redshift as expressed by the narrow emission lines. 

We have carried out a search for candidate kicked 
QSOs using spectra from the Sloan Digital Sky Survey 4 
Data Release 5 (DR5). We have focused attention on 
the broad H/3 and narrow [O in] lines, but also consider 
broad Mg n where available. Note that velocity shifts of 
broad emission lines are a well studied phenomenon often 
attributed to BLR phys i cs or even orbiting bi nary black 
holes (e.g. lGaskeUlll996tlRichards et alj|200l and refer- 
ences therein). This complicates the task of identifying 
true examples of recoil. 

2. OBSERVATIONS 

The QSOs in our study were spectroscopically iden- 
tified as QSOs in the SDSS DR5, within the redshift 
range 0.1 < z < 0.81 such that H/3 and [O ill] emis- 
sion were both measurable. The lines were measured by 
mean s of a least squares fit of a Gauss-Hcrmite func- 
tion (|Pinknev et al.l 120031 ) to the line profile, together 
with a linear fit to the continuum in the vicinity of the 
line. The broad H/3 line was fit after removing an as- 
sumed narrow H/3 line with the profile and 10% of the 
flux of the [O in] A5007 line (jBaldwin et al.|[l98l . The 
broad Mg n line, and the narrow emission lines [O n] 
and [S n] were measured when accessible. We assumed 
doublet ratios of unity for Mg n and 1.2 red-to-blue for 
[O n] . Spectral fits were accepted if the line widths and 
equivalent widths had an accuracy better than 15% and 
visual inspections showed a good fit and no artifacts. 
For more details on the m easurement procedure, includ- 
ing Fe n subtraction, see iSalviander et ail ([2007) . Our 
final data set consists of 2598 objects. The redshifts of 
the peaks of the lines were calculated from the fits to the 
line profiles. The relative displacement of the broad H/3 
line with respect to the peak of [O in] was calculated as 
Av-Hp = c(zty3 - Z[ iu])/(l + Z[ m])- The displacement 
of Mg II from [O in], AvMg, is defined analogously to 

Avu/3- 

A histogram of A«h/3 and Avyig is shown in Figure 
[TJ A Gaussian fit to the distribution of Auh/3 gives a 
mean displacement of +100 km s -1 , and the FWHM of 
the distribution is 500 km s _1 . In our sample, 40 QSOs 
(1.5%) have Ai>h/3 displacements greater than 600 km s _1 
from the mean. 

The fact that there is an overall redshift may largely 
result from physical processes in the BLR. However, we 
find an average [O in] blueshift of 30 km s^ 1 relative 
to [S n] and 40 km s relative to [O II], similar to the 
results of lBoros on ( 2005). The largest blueshifts of [O III] 
relative to [O n] were found for objects with small A«h/3- 
Interestingly, the 36 objects with Avjj/3 more than 600 
km s _1 from the mean and measurable [O n] showed 
good agreement between [O ill] and [O n]. Therefore, 
use of [O III] as a reference velocity does not significantly 
affect our calculation of the high H/3 line displacements. 

4 The SDSS website is |http: / /www.sdss.org| 
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Fig. 1. — Histogram of Ai>h/3 and Avm s along with Gaussian fits 
to the data. The distribution is somewhat broader than a Gaussian, 
with small but significant numbers of outliers with |Au| >1000 
km s . 

3. COMPARISON WITH NUMERICAL SIMULATIONS 

iSchnittman fe Buonannol (|2007l ) compute the black 
hole recoil velocities for a range of mass ratios, spin 
magnitudes and directions from post-Newtonian equa- 
tions of motion calibrated to the results of the numer- 
ical simulations. For mergers of two black holes with 
equal spin parameter a* = 0.9, and a restricted set of 
mass distributions such that rrn m^l (ml + to2) 2 > 0.16, 
ISchnittman fe Buonann o (2007) find a fraction /500 = 
0.31 of recoils greater than 500 km s _1 and a frac- 
tion /1000 = 0.079 of recoils g reater than 1000 km s" 1 . 
Convo lving the probabilities of ISchnittman fe Buonannol 
(2007) with random kick inclinations to the line of sight, 
we find the predicted observational kick fractions to be 
/ 500 = 0.18 and /iooo = 0.054. 

In our data set of 2598 objects, the fractions of large 
Av H /3 are / 50 o = 0.04 and /1000 = 0.0035. Note that 
these fractions are taken with respect to the data set as a 
whole, that is to say, assuming that every shift measured 
is due to a kick. 

4. DISCUSSION 

The observed incidence of H/3 shifts over 1000 km/s is 
several times less than theoretically expected for rapidly 
spinning holes with random orientations and similar 
masses. Moreover, closer inspection of the objects with 
large Ai>h/3 suggests that these large displacements most 
likely result from BLR physics rather than recoils. 

1. The largest shifts occur only for objects with 
large H/3 FWHM. (This di ffers from the find- 
i ngs o f ISulentic et al.1 (|2007l ) and iRichards et al.1 
(2002) for C IV A1550.) The distribution of 
FWHM for our measured H/3 lines peaks at about 
3000 km s , but for objects with shifts greater 
than 600 km s _1 from the mean we find an aver- 
age FWHM of ^5000 km s _1 , with only one ob- 
ject (SDSS J141959.21+610143.6) narrower than 
3000 km s . We have no reason to expect recoils 
to prefer larger FWHM. 

2. For recoils, all broad lines should have approxi- 
mately the same velocity shift. Figure [2] shows that 
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Fig. 2. — A plot of Mg II - and H/3 - [O III] displacements. The 
dotted line shows a fit to the data of Avyi g = 0.6 X Ai>h/3, an d 
the solid line shows the relation of equality. Removal of the high 
redshift outlier does not change the fit. 

typically At>M g ~ O.6Ai;h0 for high as well as low 
shifts. No subset of objects stands out to the eye 
as noteworthy candidates to be recoils. 

3. We find more redshifts than blueshifts, contrary to 
expectation for random kick directions. Relative 
to our measured mean, there are 59 redshifts over 
500 km s _1 and 47 blueshifts under -500 km s _1 
(out of 2598 objects). We note that this difference 
is not highly significant, and is dependent on the 
average blueshift of [O in] with respect to systemic 
velocity. 

4. We have constructed composite spectra for objects 
with Avh/3 greater than 600 km s _1 from the mean, 
less than -600 km s" 1 from the mean, and with ab- 
solute value of the shift less than 600 km s _1 In 
all three composite spectra, shown in Figure O the 
[Ne v] FWHM is several hundred km s _1 greater 
than that of the lower ionization lines. This is typ- 
ical for QSOs, where high ionization or high crit- 
ical density lines are thought to originate close to 
the black hole. The close in gas will feel the black 
hole potential and have a higher velocity co mpo- 
nent leading to a broader profile ([Laorl [20071 and 
references therein) . However, in the case of a kicked 
disk, one might expect a 'left behind' NLR or in- 
terstellar medium ionized by the off-center QSO, 
which would show line widths consistent among the 
narrow lines. Alternatively, if there were a com- 
pact NLR bound to the black hole, it would con- 
tribute a [Ne v] component centered on the velocity 
of the black hole. Neither of these cases appear 
in the composite spectra, suggesting that high- 
Auh/3 QSOs have normal NLRs. Table [T] shows 
the FWHM and equivalent width (EW) of [O n], 
[Ne in], and [Ne v] in the composite spectra of 
Figure [3j It can be seen that the shifted QSOs 
have similar narrow line intensities as the unshiftcd 
QSOs. 
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Fig. 3. — Composite spectra for the highest H/3 redshifts relative 
to [O III] (top panel, 33 objects), highest blueshifts (lower panel, 
24 objects), and shifts centering around the mean (middle panel, 
2541 objects). 



TABLE 1 

Composite spectra 



EW (A) 


redshift 


blueshift 


Av H /3l < 600 km s" 1 


[O n] 


3.07 


2.44 


2.76 


[O in] 


17.64 


19.28 


15.41 


[Ne III] 


2.73 


2.10 


1.90 


[No v] 


1.95 


1.58 


1.5 


FWHM (km s" 1 ) 


redshift 


blueshift 


Av H /3| < 600 km s" 1 


[O n] 


582 


472 


506 


[O in] 


541 


536 


420 


[Ne III] 


748 


595 


516 


[No v] 


775 


615 


650 



5. Finally, we visually inspected the high shift 
objects to see if any of them stood out as 
kick candidates. Of particular interest was 
SDSS J091833.82+315621.1, shown in Figure H 
which has the the largest H/3 shift in our sample 
= 2667 km s" 1 and Av Mg = 1231km s _1 ). 
This is very near the maximum recoil velocity 
recorded by numerical simulatio ns ([Gonzalez et al.l 
127)071 : iTichv fc Marronettill2T)07l ) . While the shifted 
H/3 line is striking in appearance and symmetric 
in shape, there seems to be little else in the spec- 
trum to distinguish it from a non-kicked BLR. NLR 
line ratios and intensities are normal, and its AvMg 
is consistent with 0.6Avnp as shown in Figure [2j 
While we cannot definitely rule out a kick, the 
shifted H/3 line in this object may just as well be 
due to some other physical process in the BLR. 

Some highly blueshifted objects, such as 
SDSS J120354.76+371137.2 (Aw H /3=-513 km s" 1 ) 
or SDSS J135800.40+404358.1 (At; H/3 =-813 
km s _1 ) have a sharp cutoff on the red side of the 
line, often accompanied by a ledge or shoulder. In 
these cases, the asymmetry shifts the peak, but 
the wings more nearly cen t er on the [O nrjredshift. 
Similarly, Richard s et al.l (|2002t ) found that large 
blueshifts (relative to [O in]) in the broad C iv 
and Mg n lines can be attributed to a diminution 
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Fig. 4.— SDSS J091833.82+315621.1, the largest shifted object 
in our sample with = 2667 km s _1 . 



of flux in the red wing of the line rather than to a 
true shift. Therefore, it seems likely that in these 
objects, the shape of broad H/3 is altered in such a 
way as to give a blueshifted peak. 

Possibly more promising are objects which, be- 
sides having symmetrical broad lines, also have 
similar Ai>h/3 and AvMg in addition to a [Ne v] 
line with similar broadening as the lower ioniza- 
tion narrow lines. Two such objects in our sam- 
ple, albeit with fairly low signal/noise spectra, are 
SDSS J134812.36+052402.6, with (Av H0 ,Avm 9 ) = 
(-706,-769) km s" 1 and FWHM ([Ne v],[0 in]) = 
(380,308) km s" 1 , and SDSS J103144.53+415420.8, 



with (Av Hp ,Av Mg ) = (-518,-462) km s" 1 and 
FWHM ([Ne v],[0 m]) = (561,572) km s" 1 . Such 
objects may merit further investigation. 



In summary, we find a number of QSOs with displaced 
broad line peaks relative to the narrow lines. However, 
for a variety of reasons, few if any of these are likely 
candidates for recoiling black holes. The retained disk 
mass in a recoil is enough to power a QSO episode last- 
ing a substantial fraction of the total QSO phase, so that 
kicked QSOs stand a fair chance of being observed if they 
occur. It is therefore likely that some mechanism is at 
work to prefer small kic ks over large ones. For example, 
iBogdanovic et all (|2007l ) propose that in gas-rich merg- 
ers, the spin and orbital angular momenta of the black 
holes become aligned with that of the large scale gas flow. 
If this process can occur before the final merger, the black 
holes will be in a configuration leading to small (<200 
km s^ 1 ) kicks, which would be undetectable through dis- 
placed broad lines. Alternatively, if black holes do not 
merge until after the QSO phase, there would be no disk 
to fuel a kicked QSO. Such a phenomenon would have to 
be detected through other means, such as tidal disr up- 
tion of stars in the galactic bulge dGezari "eTaII[200^ or 
enlargement of the galactic core ()Merritt et al.ll2004T ). 
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